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Abstract 

We discuss the sensitivity of the high-scale SUSY at 10-1000 TeV in 5°, Bg, 
K’-’ and D meson systems together with the neutron EDM and the mercury EDM. 
In order to estimate the contribution of the squark flavor mixing to these FCNCs, 
we calculate the squark mass spectrum, which is consistent with the recent Higgs 
discovery. The SUSY contribution in ck could be large, around 40% in the region 
of the SUSY scale 10-100 TeV. The neutron EDM and the mercury EDM are also 
sensitive to the SUSY contribution induced by the gluino-squark interaction. The 
predicted EDMs are roughly proportional to If the SUSY contribution is 

the level of (9(10%) for ck, the neutron EDM is expected to be discovered in the 
region of 10“^^-10“^®ecm. The mercury EDM also gives a strong constraint for 
the gluino-squark interaction. The SUSY contribution of AM/) is also discussed. 
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1 Introduction 


The supersymmetry (SUSY) is one of the most attractive theories beyond the standard model 
(SM). Therefore, the SUSY has been expected to be observed at the LHC experiments. 
However, no signals of the SUSY have been discovered yet. The present searches for the 
SUSY particles give us important constraints for the SUSY. Since the lower bounds of the 
superparticle masses increase gradually, the squark and the gluino masses are supposed to 
be at the higher scale than 1 TeV [I112113]. On the other hand, the SUSY model has been 
seriously constrained by the Higgs discovery, in which the Higgs mass is 125 GeV |1]. Based 
on this theoretical and experimental situations, we consider the high-scale SUSY models, 
which have been widely discussed with a lot of attention [5]-[20]. 

If the squark and slepton masses are at the high-scale (T(IO-IOOO) TeV, the lightest Higgs 
mass can be pushed up to 125 GeV, whereas SUSY particles are out of the reach of the LHG 
experiment. Therefore, the indirect search of the SUSY particles becomes important in the 
low energy flavor physics [2111221123] • 

The flavor physics is also on the new stage in the light of LHGb data. The LHGb 
collaboration has reported new data of the GP violation of the Bg meson and the branching 
ratios of rare Bg decays [23]-[36]. For many years the GP violation in the K and B^ mesons 
has been successfully understood within the framework of the standard model (SM), so called 
Kobayashi-Maskawa (KM) model [57], where the source of the GP violation is the KM phase 
in the quark sector with three families. However, the new physics has been expected to be 
indirectly discovered in the precise data of B^ and Bg meson decays at the LHGb experiment 
and the further coming experiment, Belle-H. 

There are new sources of the GP violation if the SM is extended to the SUSY models. The 
soft squark mass matrices contain the GP violating phases, which contribute to the flavor 
changing neutral current (FGNG) with the GP violation [38]. Therefore, we can expect the 
SUSY effect in the GP violating phenomena. However, the clear deviation from the SM 
prediction has not been observed yet in the LHGb experiment [23]-[36]. Actually, we have 
found that the GP violation of B^ and Bg meson systems are suppressed if the SUSY scale 
is above 10 TeV [32]. On the other hand, the GKMhtter group presented the current limits 
on new physics contributions of (P(10%) in i?°, Bg and systems [3D]- They have also 
estimated the sensitivity to new physics in B^ and Bg mixing achievable with 50ab“^ of 
Belle-H and 50fb“^ of LHGb data. Therefore, we should carefully study the sensitivity of the 
high-scale SUSY to the hadronic FGNG. 

In this work, we discuss the high-scale SUSY contribution to the Bg and meson 
systems. Furthermore, we also discuss the sensitivity to the D meson and the electric dipole 
moment (EDM) of the neutron and the mercury. For these modes, the most important 
process of the SUSY contribution is the gluino-squark mediated flavor changing process [31]- 
[56] . The GP violation of K meson, eK, provides a severe constraint to the gluino-squark 
mediated FGNG [571 EH]- In addition, the recent work have found that the chromo-electric 
dipole moment (cEDM) is sensitive to the high-scale SUSY [3D]. It is noted that the upper- 
bound of the neutron EDM (nEDM) [60] gives a severe constraint for the gluino-squark 
interaction through the cEDM [6T]-[66]. It is also remarked that the upper bound of the 
mercury EDM (HgEDM) [67] can give an important constraint [68] . 
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In order to estimate the gluino-squark mediated FCNC of the K, B^, Bg and D mesons, 
we work in the basis of the squark mass eigenstate with the non-minimal sqnark (slepton) 
flavor mixing. There are three reasons why the SUSY contribntion to the FCNC considerably 
depends on the sqnark mass spectrum. The hrst one is that the GIM mechanism works in the 
squark flavor mixing, and the second one is that the loop functions depend on the mass ratio 
of squark and gluino. The last one is that we need the mixing angle between the left-handed 
sbottom and right-handed sbottom, which dominates the Ai? = 1 decay processes. Therefore, 
we discuss the squark mass spectrum, which is consistent with the recent Higgs discovery. 
Taking the universal soft parameters at the SUSY breaking scale, we obtain the squark mass 
spectrum at the matching scale where the SM emerges, by using the Renormalization Group 
Equations (RGEs) of the soft masses. On the other hand, the 6x6 mixing matrix between 
squarks and quarks is taken to be free at the low energy. 

In section 2, we discuss the squark and gluino mass spectrum and the squark mixing. In 
section 3, we present the formulation of the FCNC with AF = 2 in iC, i?°, Bg and D meson 
systems together with uEDM and HgEDM. We present numerical results and discussions in 
section 4. Section 5 is devoted to the summary. The relevant formulations are presented in 
Appendices A, B, C and D. 


2 SUSY Spectrum and Squark mixing 


The low energy FCNCs depend signihcantly on the spectrum of the SUSY particles, which 
depend on the model. As well known, the lightest Higgs mass can be pushed up to 125 GeV if 
the squark masses are expected to be (T(10) TeV. Therefore, let us consider the heavy SUSY 
particle mass spectrum in the framework of the minimal supersymmetric standard model 
(MSSM), which is consistent with the observed Higgs mass. The discussion how to obtain 
the SUSY spectrum have been given in Refs. [SSUZO]. 

We outline how to obtain the SUSY spectrum in our work. The details are presented 
in Appendix A. At the SUSY breaking scale A, we write the quadratic terms in the MSSM 
potential as 

U 2 = + ml\H 2 \^ + ml{H^ • + h.c.) . (1) 

Then, the Higgs mass parameter m? is expressed in terms of m^, m\ and tan/3 as: 


m 


2 


m\ — m 2 tan^ (3 
tan^ /3 — 1 


( 2 ) 


After running down to the Qo scale, in which the SM emerges, by the one-loop SUSY Renor¬ 
malization Group Equations (RGEs) [7T], the scalar potential is the SM one as follows: 


VsM — —m^\H'^ + ttI-^ 


(3) 


Here, the Higgs coupling A is given in terms of the SUSY parameters at the leading order as 


A(Qo) = + 9 '^) 2/3 + 


1 -^ 

12 


At(Qo)-h(Qo)cot/3 
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and ht is the top Yukawa coupling of the SM. The parameters m 2 and A run with the two-loop 
SM RGBs with MS scheme [72] down to the electroweak scale Qew = and then give 

m\ = 2m?{mH) = X{mH)v^ ■ (5) 

When mn = 125 GeV is put, X{Qo) and m^(Qo) are obtained. This input constrains 
the SUSY mass spectrum of the MSSM. In our work, we take the universal soft breaking 
parameters at the SUSY breaking scale A as follows: 

= mDc{A) = mi {A) = m^c(A) = ml {i = 1, 2, 3) , 

Ml (A) = A /2 (A) = M 3 (A) = mi /2 , m\{tC) = m 2 (A) = ml , 

Au{A) = Ao2/!7(A) , A£)(A) = Ao2/z)(A) , Ae{A) = Ao2/£(A) . (6) 


By inputting mn = 125 GeV and taking the heavy scalar mass ~ Qq (see Appendix A), 
we can obtain the SUSY spectrum for the hxed Qq and tan/3. The details and numerical 
results are presented in Appendix A. 

Let us consider the squark flavor mixing. As discussed above, there is no flavor mixing 
at A in the MSSM. However, in order to consider the non-minimal flavor mixing framework, 
we allow the off diagonal components of the squark mass matrices at the 10 % level, which 
leads to the flavor mixing of order 0.1. We take these flavor mixing angles as free parameters 
at low energies. Now we consider the 6 x 6 squark mass matrix Mq in the super-GKM basis. 
In order to move the mass eigenstate basis of squark masses, we should diagonalize the mass 
matrix by rotation matrix rg' as 

m| = rg>M|rg>' , (7) 

where T^?^ is the 6 x 6 unitary matrix, and we decompose it into the 3x6 matrices as 
Tg^ = (Tgl^, in the following expressions: 


GL — 


^ GR — 



Cl3 0 

sfge *'^i3C6i 

0 

0 - 

■sf'gC 

^23 
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0 

0 -s 

23cf'3e“®'^23sge*‘^ 



^13^23^8 

0 

0 
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„R 
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^23 
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cfgcIjSee-*'^ 

nR f~tR 

'^13^23^ 

— 


cf^cgce 2 


( 8 ) 


where we use abbreviations = cos 6 *,^’^, = sinR^’^, ce = cos 9 and sg = sin 6 *. Here 

6 is the left-right mixing angle between bi and 6 /j, which is discussed in Appendix A. It is 
remarked that we take = 0 due to the degenerate squark masses of the hrst and second 
families as discussed in Appendix A. 

The gluino-squark-quark interaction is given as 


£».(«) = -iV2g.^ 5; (r“)G“ 

{<?} 






(9) 
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where L = (1 — 75 )/ 2 , R = (1+ 75 ) 72 , and denotes the gluino held, g* are three left-handed 
(1=1,2,3) and three right-handed quarks (1=4,5, 6 ). This interaction leads to the gluino-squark 
mediated havor changing processes with AF = 2 and AF = 1 through the box and penguin 
diagrams. 

The chargino (neutrahno)-squark-quark interaction can be also discussed in the similar 
way. 


3 FCNC of AF = 2 


In our previous work [39], we have probed the high-scale SUSY, which is at 10-50 TeV scale, 
in the CP violations of K, and Bg mesons. It is found that ex is most sensitive to the 
SUSY even if the SUSY scale is at 50 TeV. The SUSY contributions for the time dependent 
CP asymmetries of B^ and Bg with AB = 1 are suppresses at the SUSY scale of 10 TeV. 
Furthermore, the SUSY contribution for the b ^ sj process is also suppressed since the 
left-right mixing angle, which induces the chiral enhancement, is very small as discussed in 
Appendix A. Therefore, we discuss the neutral meson mixing P+pO(P‘^ = K, B^, Bg, D), 
which are FCNCs with AF = 2. 

In those FCNCs, the dominant SUSY contribution is given through the gluino-squark 
interaction. Then, the dispersive part of meson mixing are written as 




'12 


( 10 ) 


where are given by the squark mixing parameters in Eq.(]8|) and its explicit formu¬ 

lation is given in Appendices B and C. 

At hrst, we discuss the AB = 2 process, that is, the mass differences AMgo and AMb^, 
and the CP-violating phases 0^ and (pg. In general, the contribution of the new physics (NP) 
to the dispersive part M 72 is parameterized as 


= K2 


, (g = P°, Bg) 


( 11 ) 


where are the NP contributions. The generic hts for B^ and Bg mixing have given the 

constraints on (hg, aq) [ID], where it is assumed that the NP does not signihcantly affect the 
SM tree-level charged-current interaction, that is, the absorptive part P ^2 is dominated by 
the decay h ccs. At present, the NP contribution hg are 10-35% and 15-25% depending on 
CTg for P° and Bg, respectively. Thus, we can expect the sizable NP contribution of 0(20%). 
We will discuss whether the high-scale SUSY can hll in the magnitude of the present NP 
contribution of 0(20%). 

Next, we discuss AS = 2 process, AM^o and the CP-violating parameter in the K 
meson, ex- By the similar parametrization in Eq. fllip . the allowed region of {hx,(Tx) has 
been estimated in Ref. [40] . The NP contribution is at least 50% although there is the strong 
(Tx dependence. Therefore, it is important to examine carefully the CP violating parameter 
exi which is given as follows: 






sin 


Am(Mi^) 

V ^Mx 



^ ImA^ 


= tan ^ 


(2AMx\ 

V Ar^ ; ’ 


( 12 ) 
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with Aq being the isospin zero amplitnde in K ^ tttt decays. Here, is the dispersive 
part of the K^-K^ mixing, and AM^ is the mass difference in the nentral K meson. The 
effects of 7 ^ 0 and (j)^ < 7 r /4 give snppression effect in ex, and it is parameterized as and 
estimated by Bnras and Gnadagnoli ca as: 

= 0.92 ±0.02 . (13) 

In the SM, the dispersive part is given as follows, 

M]! = {K\n^p=2\K) 

= -^ (^) ^wBrFIMk {r]cc><lE{xc) + r]u\tE{xt) + 2r]^\c\tE{xc,Xt)) , (14) 

where Ac = VcsV*^, Xt = VtsV^*^. The E{xys are the one-loop fnnctions [71] and r]cc,tt,ct are 
the QCD corrections [73]. Then, |e^| is given in terms of the Wolfenstein parameters A, p 
and p as follows: 

|e^“| = K^C^BKlVctyX'^f] (|14bp(l - p)pttE{xt) - PccE{xc) + PctE{xc,Xt)) , (15) 


with 


GlEj^rriKM^ 

QV^tt^AMk 


(16) 


Note that |ef^| depends on the non-perturbative parameter Bk in Eq. ffTSj) . Recently, the 
error of this parameter shrank dramatically in the lattice calculations [75] . In our calculation 
we use the updated value by the Flavor Lattice Averaging Group ITS]: 


Bk = 0.766 ±0.010 


Let us write down ex as: 


, _ ,SM , 


,SUSY 
K ) 


(17) 


(18) 


where is induced by the imaginary part of the gluino-squark box diagram, which is 

presented in Appendices B and G. Since vanishes in our scheme, is given in the 

second order of the squark mixing x ^ 23 '^^ 

In addition to the above FGNG processes, the neutron FDM, arises through the cFDM 
of the quarks, dg due to the gluino-squark mixing [61]-[66]. By using the QGD sum rules, dn 
is given as 

dn = {0.79dd - 0.20du) + e{0.3dy + 0.59d^) . (19) 

where dq and denote the EDM and cEDM of quarks dq dehned in Appendix D. On the 
other hand, by using the chiral perturbation theory 

dn = e(3.0d^ ± 2.5d^ ± 0.5df) . (20) 
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Therefore, the experimental upper bound [60] 


\dn\ < 0.29 X 10-2^ecm , (21) 

provides us a strong constraint to the gluino-squark mixing. 

The HgEDM can also probe the gluino-squark mixing [ 68 ]. The QCD sum rule approach 
gives [77] 

dHg = e{d^ -d^ + 0.012df) x 3.2 x lO'^ , (22) 

and the chiral Lagrangian method gives m 

dHg = e{d^ -d^ + O.OOSldf) x 8.7 x 10"^ . (23) 

The experimental upper bound [67] 

\dHg\ < 3.1 X 10“^®ecm , (24) 

constrains the gluino-squark mixing. 

At the last step, we discuss the charm sector, which is a promising held to probe for the 
new physics beyond the SM. The mixing is now well established US] as follows: 

xd = = (3.6 ± 1.6) X 10-=" , vd = ^ = (6.1 ± 0.7) x lO’^" , (25) 

where AMd and AT o are the differences of the masses and the decay widths between the 
mass eigenstates of the D meson, respectively, and T^:) is the averaged decay width of the D 
meson. Since the SM prediction of AMd at the short distance is much suppressed compared 
with the experimental value due to the bottom quark loop, the SUSY contribution may be 
enhanced. 


4 Results and Discussions 

Let us estimate the SUSY contribution of the low energy FCNC. We calculate the SUSY mass 
spectrum at Qo = 10, 50,100,1000 TeV and interpolate the each mass of the SUSY particle in 
the region of Qq = 10-1000 TeV. This approximation is satished within 0(10%). Therefore, 
our numerical results should be taken with the ambiguity of 0(10%). The mass spectrum 
at Qo = 10 TeV is presented in Appendix A. See Refs. [39], [58] for the mass spectrum at 
Qo = 50 TeV. 

Then, we have four mixing angles 6 *^ 3 ^^ and 6 ^ 23 five phase 4>- We reduce 

the number of parameters by taking sin^h = sin^h = Sij for simplicity. In the numerical 
calculations, we scan the phases of Eq. ([H]) in the region of 0 ~ 27r for fixed s^, where the 
Cabibbo angle 0.22 and the large angle 0.5 are taken as the typical mixing. Other relevant 
input parameters such as quark masses me, m^, the CKM parameters 14^, Wb, P, V and /s, 
fx, etc. have been presented in our previous paper Ref. [55], which are referred from UTfit 
Collaboration [80] and PDG [60] . 
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Figure 1; The SUSY components of (a) AM^o and (b) AMb^ versus ttIq for S 13 = S 23 = 0.22 
(cyan) and 0.5 (blue). The horizontal red line denotes the experimental central value. 

4.1 and Bg meson systems 

At hrst, we examine the SUSY contribution in the AB = 2 process. We show the SUSY scale 
rriQ = Qo dependence of the SUSY contributions of AM bo and AMb^ in Figure 1(a) and (b), 
where the experimental central value is shown by the red line. The experimental error-bars 
are 1% and 0.1% levels for AM^o and AMb^, respectively. We take S 13 = S 23 = 0.22,0.5. 
There is no phase dependence in our predictions. It is found that the SUSY contributions 
in AMbo and AMb^ are at most 1.5% and 0.1% at ttIq = 10 TeV, respectively. Namely, the 
high-scale SUSY cannot explain the NP contributions of hd = 0.1-0.35 and hg = 0.15-0.25, 
which have been discussed in Eg. (ITT]) . As uiq increases, the SUSY contributions of both 
AMbo and AMb^ decrease approximately with the power of l/m^. Thus, there is no hope 
to observe the SUSY contribution in the AB = 2 process for the high-scale SUSY. It should 
be noted that the SM predictions are comparable to these experimental data. 

The related phenomena are the CP violation of the non-leptonic decays B^ J/ipKs 
and Bg —)■ J/^|J<p. The recent experimental data of these phases are |2Zl EH ESI ES] 

sin (f)d = 0.679 ± 0.020 , 0, = 0.07 ± 0.09 ± 0.01 , (26) 

in which the contribution of the gluino-squark-quark interaction may be included. The NP 
contributions in (f>d and (pg are expressed in terms of the parameters of Eq. flTT]) as |55j : 

(pd = 2l3d + arg(l -h , (pg = -2f3g + arg(l -h , (27) 

where (3dWs) is the one angle of the unitarity triangle giving by the CKM matrix elements 
of the SM. However, hd and hg in the high-scale SUSY are much suppressed compared with 
hd = 0.1-0.35 and hg = 0.15-0.25 of Eg. (ITT]) , one cannot hnd signals of the high-scale SUSY 
in the CP violating decays B^ —)• J/ipKs and Bg J/ip<p. 

4.2 Neutral K meson system 

At the second step, we examine the neutral K meson. We show the SUSY contributions of 
AMj^o and ex versus ttIq = Qo in Figure 2(a) and (b), where the experimental central value 
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Figure 2: The SUSY components of (a) AM^o and (b) |ex| versus uIq for S 13 = S 23 = 0.22 
(cyan) and 0.5 (blue). The horizontal red line denotes the experimental central value. 


is shown by the horizontal red line. The experimental error-bars are 0.2% and 0.5% levels for 
AMj^o and ex, respectively. Since 9 ^ 2 ^ = 0, the SUSY flavor mixing arise from the second 
order of S 13 x S 23 , where S 13 = S 23 = 0.22, 0.5 are put. 

It is found in Figure 2(a) that the SUSY contribution in AM^o can be comparable to 
the experimental value in the case of S 13 = S 23 = 0.5 whereas it is suppressed in the case 
of Si 3 = S 23 = 0.22 at rrtQ = 10 TeV. Thus, AM^o constrains the squark mixing of S 13 and 
S 23 around ttIq = 10 TeV. When the SUSY scale increases to more than 20 TeV, no SUSY 
contribution is expected. 

On the other hand, ex is very sensitive to the SUSY contribution up to lOOTeV as seen 
in Figure 2(b). The plot is scattered due to the random phases of the squark mixing. The 
experimental data of eK constrains the squark mixing and phases considerably. Actually, 
we have already pointed out that the SUSY contribution in ex could be 40% and 35% at 
ttIq = 10,50 TeV, respectively [31]. It is found that this seizable SUSY contribution still 
exist up to 100 TeV in this work. 

In the SM, there is only one CP violating phase. Therefore, the observed value of (j)d 
in Eq. fl27|) . should be correlated with ex in the SM. According to the recent experimental 
results, it is found that the consistency between the SM prediction and the experimental data 
of sin (f>d and ex is marginal. This fact was pointed out by Buras and Guadagnoli [73] and 
called as the tension between ex and sin(/)rf. Considering the effect of the SUSY contribution 
0(10%) in ex, this tension can be relaxed even if rriQ = 100 TeV. The precise determination 
of the unitarity triangle of is required in order to hnd the SUSY contribution of this level. 

It is noted that the SUSY contribution of both AMxo and ex also decrease approximately 
with the power of l/?7i^ as itiq increases up to 1000 TeV. 

4.3 The nEDM and HgEDM with ek 

The nEDM and HgEDM are also sensitive to the SUSY contribution [59l|68]. The gluino- 
squark interaction leads to the cEDM of quarks, which give the nEDM as shown in Eqs. flTOl) 
and fl2U]) . We show the predicted nEDM versus mg for the case of the QCD sum rules of 
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Figure 3: The the neutron EDM versus (a) ttIq and (b) versus for S 13 = S 23 = 0.22 

(cyan) and 0.5 (blue) for the case of the QCD sum rule. The horizontal red line denotes the 
experimental upper bound of \dn\ and the vertical one is the experimental central value of 
|ex|- 
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Figure 4: The the mercury EDM versus (a) ttIq and (b) versus for S 13 = S 23 = 0.22 

(cyan) and 0.5 (blue) for the case of the QCD sum rule. The horizontal red line denotes the 
experimental upper bound of \dHg\ and the vertical one is the experimental central value of 
\eK\- 


Eqs. flT^ in Figure 3(a), where the upper bound of |d„| is shown by the red line. The plot 
is scattered due to the random phases of the squark mixing as well as in the case of ex- 
We hnd that the contribution of EDM, dd and occupy around 25% of the neutron EDM. 
The SUSY contribution is close to the experimental upper bound up to 50TeV. Since the 
predicted nEDM depends on the phases of the squark mixing matrix signihcantly, we plot 
the nEDM versus in Figure 3(b). It is found that the predicted nEDM is roughly 

proportional to If the SUSY contribution is the level of (9(10%) for ex, the nEDM is 

expected to be discovered in the region of 10“^^-10“^®ecm. On the other hand, if the nEDM 
is not observed above 10“^®ecm, the SUSY contribution of ex is below a few %. Thus, there 
is the correlation between dn and e^^^. 
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We also show the predicted HgEDM versus rriQ for the case of the QCD sum rules of 
Eq. fl^ in Figure 4(a), where the upper bound of \dHg\ is shown by the red line. The SUSY 
contribution is close to the experimental upper bound up to 200TeV, which is much higher 
than the one of the uEDM. In Figure 4(b), we plot the HgEDM versus It is found 

that the experimental upper bound of the HgEDM excludes completely which is 

inconsistent with the experimental data. If the SUSY contribution is the level of 0(10%) for 
ex, the nEDM is expected to be discovered in the region of 10“^^-10“^®ecm. If the HgEDM 
is not observed above 10“^®ecm, the SUSY contribution of ex is below a few %. Thus, the 
mercury EDM gives more signihcant information for the gluino-squark interaction compared 
with the neutron EDM. 

However, these correlations strongly depend on the assumptions of 6^23 = and 
The deviation from these relations destroys these correlations. For instance, for the case of 
6*23 ^ ^13 with 6 *h = 9^, is much suppressed whereas the nEDM and HgEDM are still 

sizable. On the other hand, if 6 *h 3 > 9fj or 6 *h ^ 9^ is realized, the cEDMs are suppressed 
because they require the chirality flipping. In conclusion, the careful studies of the mixing 
angle relations are required to test the correlations between EDMs and 

We should comment on the hadronic model dependence of our numerical result. For both 
nEDM and HgEDM, we show the numerical result by using the hadronic model of the QCD 
sum rules in Eqs. dT^ and ([22]). We have also calculated the EDMs by using the hadronic 
model of the chiral perturbation theory in Eqs. fl2UI) and 0231) . For the neutron EDM, the 
prediction of the chiral perturbation theory is larger than the one of the QCD sum rule at 
most of factor two. However, for the mercury EDM, the prediction of the QCD sum rule is 
more than three times larger compared with the one of the chiral perturbation theory. Thus, 
predicted EDMs have the ambiguity with the factor 2 — 3 from the hadronic model. 

4.4 D-D mixing 

Since the SM prediction of AMd at the short distance is 0(10“^®) GeV, which is very small 
compared with the experimental value due to the bottom quark loop, it is important to 
estimate the SUSY contribution of AMjg. The mixing angle also appears in the up- 

type squark mixing matrix whereas the down-type squark mixing matrix contributes to 

and Bg meson systems induced by the gluino-squark-quark interaction. 

We show the SUSY component of AMd and xd versus rriq for S 13 = S 23 = 0.22,0.5 in 
Figure 5. At the SUSY scale of 10 TeV, the SUSY component may be comparable to the 
observed value. Although the accurate estimate of the long-distance effect is difficult, Cheng 
and Chiang estimated of order 10 “^ from the two body hadronic modes [SI]. This obtained 
value is consistent with the experimental one. Therefore, we should take into account the 
long-distance effect properly in order to constrain the SUSY contribution from AMd. 

Before closing the presentation of the numerical results, we add a comment on the other 
gaugino contribution. There are additional contributions to the FCNC induced by chargino 
exchanging diagrams. The chargino contribution to the gluino one is approximately 10% in 
the above numerical study of AF = 2. Thus, the chargino contributions are the sub-leading 
ones. 
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Figure 5: The SUSY component of {&)AMd and (b)x£) versus ttIq for S 13 = S 23 = 0.22 (cyan) 
and 0.5 (blue). The horizontal red line denotes the experimental central value. 


5 Summary 

We discussed the sensitivity of the high-scale SUSY at 10-1000 TeV in the 5°, Bg and 
meson systems. Furthermore, we have also discussed the sensitivity to the D-D mixing, the 
neutron EDM and the mercury EDM. In order to estimate the contribution of the squark 
flavor mixing to these FCNC, we calculate the squark mass spectrum, which is consistent 
with the recent Higgs discovery. 

The SUSY contributions in AM^o and AMq^ are at most 1.5% and 0.1% at tUq = 10 TeV, 
respectively. As ttIq increases, the SUSY contributions of both AM^o and AMb^ decrease 
approximately with the power of 1/m?^. Therefore, the SUSY scale increases to more than 
10 TeV, no signal of the SUSY is expected. On the other hand, the SUSY contribution in 
AMxo can be comparable to the experimental value in the case of S 13 = S 23 = 0.5 whereas 
it is suppressed in the case of S 13 = S 23 = 0.22 at ttiq = 10 TeV. Furthermore, the SUSY 
contribution in ex could be large, around 40% in the region of the SUSY scale 10-100 TeV. 
By considering the effect of the SUSY contribution 0(10%) in ex, the tension between ex 
and sin0d can be relaxed even if the SUSY scale is 100 TeV. 

The neutron EDM and the mercury EDM are also sensitive to the SUSY contribution 
induced by the gluino-squark interaction. The \dn\ is expected to be close to the experimental 
upper bound even if the SUSY scale is 50 TeV. The predicted nEDM is roughly proportional 
to If the SUSY contribution is the level of 0(10%) for ex, the \dn\ is expected 

to be discovered in the region of 10“^^-10“^®cm. For the {dxgl, the SUSY contribution is 
close to the experimental upper bound up to 200TeV, which is much higher than the one of 
the nEDM. If the HgEDM is not observed above 10“^®cm, the SUSY contribution of ex is 
below a few %. Thus, the mercury EDM gives more signihcant information for the gluino- 
squark interaction compared with the neutron EDM. It may be important to give a comment 
that these predictions depend strongly on the assumptions of ^23 = ^13 ^J^cl The 

deviation from these relations destroys these correlations. In conclusion, the careful studies 
of the mixing angle relations are required to test the correlations between EDMs and 
The predicted EDMs have also the ambiguity with the factor 2 — 3 from the hadronic model. 
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Since the SM prediction of AMo at the short distance is 0(10 GeV, which is very small 
compared with the experimental valne, it is important to estimate the SUSY contribntion of 
AMd. 

In conclusion, the more detailed studies of meson system, the EDMs of the neutron 
and mercury are required in order to probe the high-scale SUSY at 10-1000 TeV. 
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Appendix A : Running of SUSY particle masses 

In the framework of the MSSM, one obtains the SUSY particle spectrum which is consistent 
with the observed Higgs mass. The numerical analyses have been given in Refs. [691110]. At 
the SUSY breaking scale A, the quadratic terms in the MSSM potential is given as 


U 2 = + ml\H 2 \^ + ml{Hi ■ H 2 + h.c.) . 

The mass eigenvalues at the Hi and H 2 = ^H^ system are given 


(28) 


9 m? + rrin 

mi- = — - - - =F 

T r) ^ 


mf — m 


2 \ I 4 


(29) 


Suppose that the MSSM matches with the SM at the SUSY mass scale Qq = mo. Then, the 
smaller one mii_ is identihed to be the mass squared of the SM Higgs H with the tachyonic 
mass. The larger one m^ is the mass squared of the orthogonal combination H, which is 
decoupled from the SM at Qo, that is, m^ ~ Qq. Therefore, we have 


m^ = —m‘^{Qo) , m^ = m^(Qo) = m^ + m^ + m^ , 


with 


4 /2| 2\/2i 2\ 

m 3 = [mi + m )[m 2 + m ) , 

which leads to the mixing angle between Hi and H 2 , P as follows: 


(30) 


(31) 


9 ^ 

tan /3 = —9 -^ 

mi + m^ 


H = cos l3Hi + sin f3H2 , "H = — sin (]Hi + cos I 3 H 2 ■ (32) 


Thus, the Higgs mass parameter m? is expressed in terms of m\, m\ and tau/S: 

2 ml-mliw?l 3 

m = -;— . (33) 

tan^ (3-1 ^ ^ 

Below the Qq scale, in which the SM emerges, the scalar potential is the SM one as follows: 


VSM ^ -iAh? + ^\H 


(34) 


Here, the Higgs coupling A is given in terms of the SUSY parameters at the leading order as 


KQa) = -Ag + 9 ' )cos 2/3 + 


87r2' 


1 - — 
12 


At[Q q) - ^i[QQ) coi (3 

X,- --. (35) 


and ht is the top Yukawa coupling of the SM. The parameters m 2 and A run with the SM 
Renormalization Group Equation down to the electroweak scale Qew = ^Hi and then give 


m^j = 2m^(m//) = X[mH)v^ 


(36) 
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It is easily seen that the VEV of Higgs, (H) is v, and {H) = 0, taking acconnt of (Hi) = v cos (3 
and {H 2 ) = V sin where v = 246 GeV. 

Let us £x rrifj = 125 GeV, which gives A(Qo) and m^(Qo)- This experimental input 
constrains the SUSY mass spectrum of the MSSM. We consider the some universal soft 
breaking parameters at the SUSY breaking scale A as follows: 

= m£,o{A) = = m^c(A) = ml {i = 1, 2, 3) , 

Mi(A) = M 2 (A) = M3(A) = mi /2 , ml{A) = ml{A) = ml , 

= Aoyu{A) , A£)(A) = Ao|/d(A) , Ae{A) = Ao 2 /£(A) . (37) 

Therefore, there is no flavor mixing at A in the MSSM. However, in order to consider the 
non-minimal flavor mixing framework, we allow the off diagonal components of the squark 
mass matrices at the 10% level, which leads to the flavor mixing of order 0.1. We take these 
flavor mixing angles as free parameters at low energies. 

Now, we have the SUSY hve parameters. A, tau/S, mg, mi/ 2 , Aq, where Qq = mg. In 
addition to these parameters, we take fi = Qq. Inputing mn = 125 GeV and taking ~ Qq, 
we can obtain the SUSY spectrum for the hxed Qq and tau/S. 

We present the SUSY mass spectrum at Qq = 10 TeV. The input parameter set and 
the obtained SUSY mass spectra at Qo are summarized in Table 1, where we use rnt{mt) = 
163.5 ± 2 GeV [ 6 OI l80] . These parameter sets are easily found from the work in Ref. [69]. 


Input at A and Qo 

Output at Qo 

at A = 10^^ GeV, 

mg = 12.8 TeV, m^ = 5.2 TeV, m^ = 2.9 TeV 

mo = 10 TeV, 

m^^ = mj^ = 12.2 TeV 

mi /2 = 6.2 TeV, 

mi^ = 14.1 TeV, m^^ = 8.4 TeV 

Ao = 25.803 TeV; 

= 15.1 TeV 

at Qo = 10 TeV, 

^sn,dR - ^CR,UR = 14.6 TeV, m^ = 13.7 TeV 

/i = 10 TeV, 

mf^ = m^^^ = 10.4 TeV, mf^ = 9.3 TeV 

tan = 10 

^RL,eL = = 10-8 TeV, = 10.3 TeV 

Xt = -0.22, Xh = 0.126 


Table 1: Input parameters at A and the obtained SUSY spectra at Qo = 10 TeV. 


As seen in Table 1, the hrst and second family squarks are degenerate in their masses, 
on the other hand, the third ones split due to the large RGB’s effect. Therefore, the mixing 
angle between the hrst and second family squarks vanishes, but the mixing angles between the 
hrst-third and the second-third family squarks are produced at the Qo scale. The left-right 
mixing angle between and bji is given as 


mb(Ab(Qo) -/itan,d) 

o' ~-5-5- 

mi — mi 

ol Ofl 


(38) 


which is very small, 0(0.01) at 10 TeV. The lightest squark is the right-handed stop and the 
lightest gaugino is the Bino. 
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Appendix B : Squark contribution in AF = 2 process 

The AF = 2 effective Lagrangian from the gluino-sbottom-quark interaction is given as [82] : 


^AF=2 

^eS 


[CvllOvll + CvrrOvrr] 


E 

i=l 


/n(h 


+ C 


(i) 

SRR^SRR 


+ c 


W n^i) 

SLR^SLR 


(39) 


where 

OvLL = {QaF^Q°‘){QbFLQ’^), OvRR = {qaFRQ^) ■, 

= {qaLQ^WbLQ^). = {qaLQ^){q,LQ^), 

O^sIr = (qaRQ’^MRQ^), = {qaRQ^){q,RQ‘^), 

4lh = (gaEQ“)(®i?Q'), = (qaLQFqbRQ^), (40) 

with (P, Q, g) = (- 8 °, &, c?), {Bs,b,s), {K^,s,d). The L, P denote (1 ± 75 ) 72 , and a, 6 are color 
indices. Then, the P^-P^ mixing, M 12 , is written as: 


Mi2 — 


1 

2mp 


(pVf/=^ 


p°). 


(41) 


The hadronic matrix elements are given in terms of the non-perturbative parameters P* as: 


{P°\Ovll\P") = '^rnlflBi, {P«\Ovrr\P°) = {P°\Ovll\P°), 

{P»|C>SlI^°) = -RrPpflRpB,, (-P“|OS,±l-P”) = (-P'IOSlI^”). 
(P»|OsllI^°) = ^mUlRpB:,, (PVsrhI^") = {^’"|oSl7“>. 

(F"\OfRR\F) = \mlSlRpBt, (P“|0?h|.P°> = \m%SlRpB,, (42) 


where 


Pp — 


mp 


mq ± niq^ 

The Wilson coefficients for the glnino contribution in Eq. 


(43) 


are written as 
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CvLli'^g) 

CvRni.'mg) 

^SRRi'^g) 

^SLli'^g) 

^SRRi^'^g) 

^slhi^'^g) 


E(a 


2 

mi 

9 IJ=l 


{d) YJ (\(d) \ij 
gllIi K'^gllIj 




— CvLL{xng){L -H- i?), 

(^GLr)I i^GLRJ J 


a 

m 


s 
2 

9 /,J=1 
^(1) 


{d) \ij(\{d) 

9 


^i[i](^?,^j), 


~ ^SRRi'^g)i^ ^ 
6 


a 

m 


2 
s 
2 

9 I,J=1 
72) 



GLr)i {^GLr}] ( 2 ) 


~ ^SRR('^g)(^ ^ 


^SLk("^<7) 


2 

mi 


GLRJI V'GRLU \ J 92[1]{x^i-,X^j) 


I V KjL/it/1 \ \jii\,L//J 

/,J=1 I 

_i_ ( \id) Y9(\id) Yi 
+ K'^GLUI K'^grrJj 


11\ 


yC/l[l](xf,x5) - ^^2[1 ](x?,x5) 




Y^ |(•^GLi^)/(-^Gi^L)i ^ 92[1]{x\,X^j) 


\(\id) Yj(\{d) \*i 


2 - - ID - - 

g^i[i](a;?,a;5) + —92iiM,x^j) 


where 


(AgL)l = (rLT)f(rgl) 


GLy'A' ; 


(\(d) \*i /p(“)t\-fS'/p(“) \ 

~ I-*- GL )i I-*- GR)K ) 


(\Kd) Yd — V 

y-^GRR/K — I-*- GRti I-*- GR)K > 

{ \id) Yd _ /p('^)t\-f4'/p('^) \ j 

K'^GRLIK — I-*- GRii I-*- GL)K ■ 


(44) 


(45) 


Here we take {i,j) = (1,3), (2,3), (1,2) which correspond to 5°, Bg, and mesons, 
respectively. The loop functions are given as follows: 


If xf ^ (x^j = m] /ml), 




Xj,Xj) = 


X^J log Xj 


92 [ 1 ]{ 


Xj,Xj) = 



Xj 


Xj log Xj 


+ 


X-^: 


[x- 


.g \2 


log X® 


xt 


X" 


1)= 


+ 


X ) 


(46) 
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2 


If xf = Xj, 


^l[l](xf,xf) = 
^2[l](xf,xf) = 


XI 


1) logxf 


+ 


(xf - 1)3 (xf - 1)2 ’ 
2 xf log xf ^ xf + 1 


(xf - 1)3 (xf - 1)^ 


(47) 


Taking account of the case that the gluino mass is much smaller than the squark mass scale 
Qo, the effective Wilson coefficients are given by using the RGEs for higher-dimensional 
operators in Eq. flS^ at the leading order of QCD as follows: 


CvLii'iTT'bi-^ — 2 GeV)) —r]y^LLCvLL{Qo), CvRR{rnb{A — 2 GeV)) — riy)^jlCvLL{Qo), 

= 2 GeV))t 
cfl,{m,(A = 2 GeV))) 

= 2 GeV))) 

= 2 GeV)) / 




^sliiQo) 

C-SLliQo 


c 


( 1 ) 


^SRR\ 


srr^Qo) 
^(2) 

'^SRR 

hi) 


{Qo 


X^lC^^XeL, 


XrR^RR ^Xrr, 


ci}Um,(A = 2 GeV))) ^ (C),‘'r(Q, 

p(2) 


45H(mi(A = 2 GeV)) ) 


'SLR 


(Qo)J 


where 


VvLL — VvRR — 


XlWef'XRR, 


as(y»i) Y 

as{mb)) 


(48) 





0-1 

^LRi 


VVLL — VvRR — 


as{mb) 


as{mt) J \as{mb) J \as{A = 2 GeV) 
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Vll — Vrr — SlL 




VAg = 



Vlr — ^lr 

'as{mt) 


'WV 24 T 16 - 
SlL = ( 60 60 


^LL — ^F/R — 


RR 



'lAq 


0 


as{mb) 


a,(A = 2 GeV), 

d^R = —16, 


^2 _ 2 

“lh — 


(49) 


For the parameters {i = 2 — 5) of B mesons, we use values in [S3] as follows: 


5f4)(^^) = 0.79(2)(4), 

Bf'tmj) =0.92(2)(4), 


Bf‘^)(m,) = 1.15(3)(15), 

= 1.72(4)(if), 


= 0.80(1)(4), 

Bf‘’(mi.) = 0.93(3)(8), 


fif“)(m,) = 1.16(2)(16), 

Bf'>(mi) = 1.75(3)(3‘) . 

(60) 

On the other hand, we use the most updated values for b[^"' and b[^^ as [80]: 


= 1.33 ± 0.06 , 


(61) 

For the paremeters B^Y = 2 — 5), we use 

1 following values [84]. 


Bf\2GeY) = 0.66 ±0.04, 

Bf\2GeY) = 1.05 ±0.12, 

(52) 

Bf\2GeY) = 1.03 ±0.06, 

Bf\2GeY) = 0.73 ±0.10, 


and we take recent value of Eq. ffT7|) for deriving i?|^^(2GeV). 

For the paremeters i?f (f = 1 — 5), we use following values [85l 186] . 

b[^^ (3GeV) = 0.75 ± 0.02, (3GeV) = 0.66 ± 0.02, (3GeV) = 0.96 ± 0.05, 
fii^^(3GeV) = 0.91 ± 0.04, 5j^^(3GeV) = 1.10 ± 0.05. (53) 


Appendix C : The loop functions Fi 


The loop functions Fi{x^) are given in terms of = m^/m^ (/ = 3, 6) as follows: 


r xi log xi 

F, (xl) = ^ ^ 

2{xl-lY 


+ 


x^gF — 5x1 — 2 
12(xl - 1)3 


I^ H)^logxl 2 (x 1)2 +5x1-1 

F2{Xg) = -7T# + 


2(xl - 1)4 


12(xl - 1)3 


FGxi) = 


logxl 


+ 


Xt, 


xi +1 


(xl - 1)3 2(xl - 1) 


xpogxl -3 , ^ t I F 

F) = -pppi + 2(F^ ^ ■ 

(54) 
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Appendix D : EDM and Chromo-EDM of quarks 

We present the EDM of the strange quark from the gluino contribution as the typical example 


MQa) = - 2 ^ 47Ta(mj)Im[A2^‘(Qa)l, 


(55) 


where 




Q^a^(mg) 8 ^ 1 


47r 


3 ^ 2m', 

1=1 dj 


+ ^s{\fRR)f) (^2^) 


(^4^) 


(56) 


On the other hand, the chromo-EDM (cEDM) of the strange quark from gluino contribution 
is given as: 

d^{.Qo) = -2^47ra,(mg)Im[klf^(go)], (57) 

where 


Af\Q«) = 


a,(mj) 1 


6 




1 


, + ^s{\GRR)f ) (9Ei(a:f) + F2(xf) 

7=1 “""di 

(d) \22 


Ad) x22 






(58) 

( 69 ) 


Including the RGE effect of QCD [87], the cEDM of the strange quark is given as 
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25 


(60) 


On the other hand, the EDM operator is mixied with the cEDM operator during RGE 
evolution. Then, one obtains 

16 16 m 

ds(2GeV) — ds{Qo) ( — j —r ) ( — j r ) ( — ... ) H- d^ (Qo) x (61) 

\as{mt)J \as{mb) J \as{2GeV) J 


as{^\ 21 / as(^A / asjmb) \ _ / asjQo \ / as(^A / Qs(^fe) 

as{mt)) \as{mb)) \as{2GeY)) \as{mt) J \as{mb) J V“?*(2GeV) 

The EDMs and cEDMs of the down- and up-quarks induced by the gluino interaction are 
also given by the similar formulas. 
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